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Abstract 

Motivated by the large discrepancy of CP-violating phase in B s — B s mixing between 
the experimental data and the Standard Model prediction, we pursue possible solutions 
within a family non-universal Z' model. Within such a specific model, we find that 
both the B s — B s mixing anomaly and the well-known "irK puzzle" could be moderated 

simultaneously with a nontrivial new weak phase, $ 72° (SI) or -82° (S2). With 

the stringently constrained Z 1 coupling B^ b , we then study the Z 1 effects on the rare 
B — > X s /j, + ij,~ and B s — > [i + fi~ decays, which are also induced by the same b — > s 
transition. The observables of B — > X s [i + n~ , at both high and low q 2 regions, are found 
to be able to put strong constraints on the // — // — Z' coupling, B^ ~ 10~ 2 . It is also 
shown that the combined constraints from B s — B s mixing, B — > irK and B — > X s [.i + fi~ 
do not allow a large Z' contribution to the pure leptonic B s — > fj- + fi" decay. 
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1 Introduction 



As particle physics is entering the era of LHC, one may expect direct evidences to be available 
to establish whether new particles and interactions are present. Meanwhile, high sensitivity 
studies of low energy phenomena would complement the direct discovery physics at LHC. The 
flavor changing neutral current (FCNC) processes, such as b — > s transitions, arise only from 
loop effects within the Standard Model (SM), and are therefore very suitable for testing the SM 
and probing its various extensions. Recently, both CDF and DO collaborations have announced 
the measurements of CP violation parameters in B s system, with the obtained CP-violating 
phase 

S = -0.57i£H(stat)i£E£(syst) DO collaboration [1] , (1) 
<p s E [-2.82, -0.32] (68%C.L.) CDF collaboration [2] , (2) 

while within the SM this phase is expected to be 

0f M = -2/3f M = -2 nrg[-V ts V t l/(V cs V: b )} = -2 x (0.018 ± 0.001) , (3) 

which deviates from the DO measurement Eq. ([1]) by more than 2a. Combining all the available 
experimental information on B s — B s mixing, the UTfit collaboration claims that the divergence 
of 4> s between the experiment measurements and the SM prediction is more than 3a [3J. Taking 
into account the deviation of <p s in a generic scenario of NP, the CKM-fitter group has found that 
the SM is disfavored at 2.5<r [I]. Interestingly, comparing an updated theoretical predication 
of B s — B s mixing with DO [S] and CDF [B] early results based on lfb' 1 data, the authors 
of Ref. [7] have found the mixing phase 2a deviated from the SM expectation. Such a large 
observed phase, if still persisting in the upcoming experimental measurements, would indicate a 
signal of new physics (NP) manifested in b — > s transitions. In the following numerical analyses, 
we would use the UTfit results of 4> s [3] as benchmarks. 

Motivated by the above observed anomaly, in this paper we shall pursue possible solutions 
within a family non-universal Z' model [8], which could be naturally derived in certain string 
constructions [9], E e models [TU] and so on. Searching for such an extra Z' boson is an important 
mission in the experimental programs of Tevatron [TT] and LHC [T2]. Performing constraints 
on the new Z 1 couplings through low-energy physics is, on the other hand, very important and 
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complementary for direct experimental searches. It is interesting to note that, within such 
a specific scenario, both the CP- violating phase problem and the well-known u ttK puzzle" in 
hadronic B — > irK decays could be resolved [T3| [HI 03] • Since both the B s — B s mixing and the 
B — > tvK decays involve the same b — s — Z' couplings, it is worthwhile to perform a constraint 
on these couplings with all the available experimental data taken into account simultaneously 
At the same time, we could also get the allowed ranges for flavor- conserving u — u — Z' and 
d — d — Z' couplings. 

The FCNC b — > sl + l~ {I = e, /i, r) transition, which gives rise to the rare inclusive B — > 
X s fi + fi~ and the purely leptonic B s — > /i + yU~ decays, is another important process to probe 
NP. Averaging the recent experimental data from BABAR |16j . Belle [17] and CLEO[18], the 
Heavy Flavor Averaging Group (HFAG) presents the following total branching ratio [TH] 

B(B -> X s fi + fi-) = (4.3±i;l) x 10" 6 . (4) 

As for the ones in the low (1.0 GeV 2 < q 2 < 6.0 GeV 2 ) and high (14.4 GeV 2 < q 2 < 25GeV 2 ) 
q 2 regions, after naively averaging the BABAR [16J and Belle [TTJ measurements, we get re- 
spectively 

B L (B -> X s n + yr) = (1.6 ± 0.5) x 10" 6 , (5) 
B H (B X s fi + /2-) = (0.44 ± 0.12) x 10" 6 . (6) 

Theoretically, with the up-to-date input parameters, the SM predictions [201 EI] for the above 
three observables are about 5.0 x 10~ 6 , 1.8 x 10 -6 and 0.45 x 10~ 6 respectively, which agree 
with the experimental data well. It implies that such observables in B — > X s fi + fi~ , together 
with the measurements of B s — B s mixing and hadronic B — > irK decays, may provide strict 
constraints on the new Z' couplings involving the lepton sector. 

As for the B s — > decay, in addition to the electro-weak loop suppression, the decay 

rate is helicity suppressed in the SM and predicted to be about 3 x 10~ 9 [201 E2J |23], which is 
still one order of magnitude lower than the CDF upper bound [21] 

B(B S -> fi + /2~) < 4.7 x 10~ 8 (90% C.L.). (7) 

It is expected that precise measurements would be available at the upcoming experiments at 
LHC and super B factories. As a consequence, we shall also investigate the Z' contribution to 
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this decay mode within the parameter spaces constrained by B s — B s mixing, B —t ttK and 
B — > X s fi + fi~ decays. 

Our paper is organized as follows. In Section 2, after a brief review of B q — B q mixing 
within the SM, we pursue possible solutions to the B s — B s mixing anomaly within a family 
non- universal Z' model, taking into account the constrains from B — > irK decays [15]. In 
Section 3, the effects of such a NP scenario on B — > X s /i + /i~ and B s — > fi + fi~ decays are 
investigated in detail. Our conclusions are summarized in Section 4. Appendix. A includes all 
of the theoretical input parameters. 

2 Constraints on Z' couplings from B q — B q mixing and 
B — > ttK decays 

2.1 Theoretical framework 

Within the SM, the effective Hamiltonian Ufff(AB = 2) for B q -B q mixing, relevant for scales 
fib = 0(rrib) is given by [20] 

Gj 

167T 2 



H s e f f (AB = 2) = -/- 2 M 2 w {V th V; q ) 2 C Q {fx b )Q{AB = 2) + h.c. , (8) 



where Q(AB = 2) = (qb)v-A(lb)v-A- Accurate to next-to-leading order (NLO) in QCD, the 
off-diagonal term M^ M (g) in the neutral B-meson mass matrix is given by 

2m B M 1 s 2 M (q) = (B q \%lf f {AB = 2)\B q ) 

= ^M 2 w (V tb V*)\B B J 2 Bq )ml qVB S (x t ) , (9) 

where Mw is the mass of W boson, B Bq the "bag" parameter, and f Bq the B-meson decay 
constant. Explicit expressions for the short- distance QCD correction function rj B and the 
"Inami-Lim" function 5*o(x t ), with x t = mt j^i^ , could be found in Ref. |2U] . 

Recently UTfit collaboration has performed a model-independent analysis of NP effects to 
B q — B q mixing in terms of two parameters C Bq and <p Bq , with the following parametrization [3] 



in \yf ull \B ) asm i4>l M , a np i(2^+^f M ) 

~ (B q \n s ji\B q ) A sM e i*r 



Table 1: Fit results for B q — B q mixing parameters Cs q and <pB q by UTfit collaboration [3J. The 
two solutions for <pB s , SI and S2, result from measurement ambiguities, see Ref. [3] for details. 



NP parameters 




M] 


c Bs 


Bs [°] (S2US1) 


68% prob. 


0.96 ±0.23 


-2.9 ± 1.9 


1.00 ±0.20 


(-68.0 ±4.8) U (-20.3 ±5.3) 


95% prob. 


[0.57, 1.50] 


[-6.7,1.0] 


[0.68,1.51] 


[-77.8,-58.2] U [-30.5,-9.9] 



Within the SM, the modulus Cs 9 and the phase 4>B q are predicted to be one and zero, respec- 
tively. Combining all the available experimental information on B q — B q mixing, the fitting 
results at 68% and 95% probabilities from Ref. [3] are listed in Table. [TJ For each probability, 
UTfit has found two solutions for (f)B s due to measurement ambiguities [3J: one is close to, but 
still 3a deviated from, the SM expectations (denoted as SI hereafter); another one is much 
more distinct from the SM and even require dominant NP contributions (S2). Such large devi- 
ates may suggest the first evidence of NP exhibited in b — > s induced processes[3] . So, in the 
following we pursue possible solutions within a family non-universal Z' model [8]. 

While the general framework for Z'-induced FCNC transitions has be formulated by Lan- 
gacker and Plumacher [8], our discussion throughout this paper for the Z' contributions goes 
under the following simplifications: (1) neglecting kinetic mixing since it only amounts to a 
redefinition of the unknown Z 1 couplings; (2) neglecting the Z — Z' mixing, which has been 
constrained to be tiny by the Z-pole measurements at LEP [25j [26], but can be easily incor- 
porated [SI EZj ; (3) no significant renormalization group (RG) evolution effects between Mz> 
and My/ scales; (4) although there are no severe constraints on right-highted q — q — Z' cou- 
plings, we follow the simplification in the literature [131 EH EH] and assume that right-handed 
couplings are flavor- diagonal and hence real due to the hermiticity of the effective Hamiltonian, 
while flavor-off- diagonal left-handed coupling terms will result in sizable FCNC — Sl — Z' 
couplings. 

Then, the effective Hamiltonian UKAAB = 2) induced by Z' contribution at M\y scale 
could be written as 

Hf ff (AB = 2) = %i^) 2 Q(AB = 2) + h.c. , (11) 
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where Bk is the Z' — b — q coupling, whose definition is different from the one used in our 
previous paper [15] by a factor 9 ^^ z , with g\ and being the gauge couplings of Z and 
Z' bosons, respectively. Due to our assumed simplifications, the RG running of the Wilson 
coefficient induced by Z' boson is the same as that of the SM, with the corresponding evolution 
matrix U LL (fi b , M w ) given to the NLO level by [20] 



U LL ( f i b ,Hw) = [1 + ^J 5 ]C/V,^)[1-^^J 5 ], (12) 

with U°(fib,^w) = i a s{^w) / a s{^b)) 2l3<) i 7q = 4, 0o — 23/3, and J 5 = 1.627 in naive dimensional 
regularization (NDR) scheme with 5 effective quark flavors. 

After some simple derivations, one can get the final Z' contribution to M\2(q) 

2m Bq M*{q) = ^^ L |^| 2 e^^m| 9 (5 B9 /i 9 ) , (13) 

with 

(0) 

U' LL ^(a s ^ w ))^[l-^lj 5 }. (14) 
Finally, we get the total contribution to the off-diagonal mass matrix term 

M 12 (q) = M? 2 M (q) + M? 2 \q). (15) 

The mass difference, which describes the strength of the B q — B q mixing, is then given by 
AM q = 2|M 12 (g)|. An early general investigation of Z 1 effects in B q — B q mixing could be 
found in Ref.[27]. 

2.2 Numerical results and discussions 

With the UTfit results at 68% and 95% probabilities listed in the Table. [T] as constraints, 
respectively, we get the allowed ranges for the Z' parameters as shown in Fig. [H with the 
corresponding numerical results given in Table. |5J We find that the new b — s — Z' coupling, 

with a new weak phase <\> L S ~ -58° (<f% 80°) and strength \B^ b \ ~ 1.2 x 10~ 3 (2.2 x 10~ 3 ) 

corresponding to the UTfit result SI (S2), is crucial to resolve the observed B s — B s mixing 
phase anomaly. On the other hand, the AM^ is well measured and in good agreement with 
the SM predictions, the strength of b — d — Z 1 coupling involved in B^ — B® mixing should be 
much weaker than the one of the SM box diagrams. Numerically, \B% b \ is found to to be about 
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Figure 1: The allowed regions for the parameters \Bh\ and 4>^ under the constraints from 
B q — B q mixing and B — > ttK decays. SI and S2 correspond to the two solutions in Table 1. 



Table 2: Numerical results for the parameters \Bk\ and 0^ under the constraints from B q — B q 
mixing only. For each solution, the first (second) row corresponds to the constraints from Cb 
and 4>B q at 68% (95%) probability. 



Solutions 


l^ 6 |(xio- 3 ) 




\BH(xlO- 3 ) 




SI 


1.24 ±0.16 


-58 ±6 


0.12 ±0.03 


-23 ±21 




1.18 ±0.29 


-52 ± 13 


^ 0.26 


arbitrary 


S2 


2.17±0.07 


-80 ±2 








2.19 ±0.14 


-80 ±4 







1.2 x 10~ 4 . Combining with the constraints by B® — B® mixing and B — > nK decays, we find 
that the relative strength, \B^ b /B^ b \ ~ O(10 _1 ), is quite similar to the hierarchy of CKM matrix 
elements within the SM, |V^Vt&/V^*Vt&| ~ 0.2. On the theoretical side, it is noted that such a 
hierarchy is not required by the Z' model itself. Although the Z' model considered here is not a 
model of Minimal Flavor Violation (MFV) type, the low energy effective Hamiltonian resulting 
from the flavor-changing Z' couplings also obey the so-called MFV hypothesis [291 ED] driven by 
the experimental data available so far, which may imply that the flavor-changing interactions 
in the Z' model are also linked to the known structure of the SM Yukawa couplings. 

From Eqs. ([9]), fllOp and (|13j) . we find that the parameters Cs q and 4>B q are independent of 
the theoretical uncertainties associated with the non-perturbative factor BB q fs q within such a 



Table 3: Numerical results for the mass difference AM q (ps 1 ). 



Solutions 


Exp. dg 


SM 


SI 


S2 


AM, 


0.508 ±0.005 


0.525 ±0.057 


0.522 ±0.077 




AM S 


17.77 ±0.12 


18.18 ±1.47 


17.14 ±2.30 


17.42 ± 2.40 



family non-universal Z' model under our assumed simplifications. However, to get the mass 
difference AM,, such uncertainties are unavoidable. With the relevant input parameters listed 
in the Appendix A, the final numerical results for AM q are listed in Tabled It can be seen that, 
after including the Z 1 contributions, our predictions for AM 9 also agree with the experiment 
data, taking into account the respective theoretical uncertainties. 

In our pervious paper [15], we found that a nontrivial new weak phase 0f ~ —86° associated 
with the b — s — Z 1 coupling is helpful to resolve the so-called u irK puzzle" , which is similar to 
our present fitting result 0f ~ —58° in SI (</>^ ~ —80° in S2) from B s — B s mixing. However, 
as found in Ref. [15], the range 0f > —50° is almost excluded by the CP-averaged branching 
ratios and direct CP asymmetries of B — > ttK decays. So, it is very necessary and interesting to 
re-evaluate the ranges of Z' couplings under the constraints from B s — B s mixing and B — > 7rK 
decays simultaneously. 

Like the Case IV in Ref. [15], we give up any simplifications on the flavor-diagonal u — u — Z' 
and d — d — Z' couplings, and use the QCD factorization (QCDF) |3TJ approach to calculate 
the amplitudes of B — > ttK, ttK* and pK decays. As for the end-point divergence appearing 
in twist-3 spectator and annihilation amplitudes, instead of the parametrization scheme, we 
quote an infrared finite dynamical gluon propagator derived by Cornwall [32] to regulate it. 
Explicitly we quote m g = 0.50 ± 0.05 GeV, which is a reasonable choice so that most of the 
observables for B — > irK, ttK* and pK decays are in good agreement with the experimental 
data [33]. In this way, we find that the time-like annihilation amplitude could contribute a 
large strong-interaction phase, while the space-like spectator-scattering amplitude is real [33J. 
The explicit comparison for the two schemes have been systemically discussed in our previous 
papers [15l [33] . Although numerically these two schemes have some differences, both of their 
predictions are consistent with most of the experimental data within errors. 
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Figure 2: The allowed regions for the parameters B^ dd under the constraints from Cb s , 
4>b s (95% prob. only) and B — > irK decays. 
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Table 4: Numerical results for the parameters \B^ b \, B uudd and 0^ under the constraints from 
B s — B s mixing and B — > irK decays. The other captions are the same as the ones in Table [2J 



Solutions 


|i?il(xl0- 3 ) 






B R 


B dd 


n dd 


SI 


1.18 ±0.16 


-62 ±5 


0.66 ±0.38 


-0.13 ±0.12 


0.88 ± 0.36 


-0.18 ±0.14 




1.09 ±0.22 


-72 ±7 


0.34 ±0.55 


-0.04 ±0.18 


0.70 ± 0.48 


-0.07 ±0.20 


S2 


2.19 ±0.06 


-81 ±2 


-0.02 ± 0.34 


0.01 ±0.12 


0.22 ±0.32 


0.01 ±0.12 




2.20 ±0.15 


-82 ±4 


0.02 ±0.34 


-0.01 ±0.12 


0.27 ±0.32 


-0.04 ±0.24 



Including the constraints from B s — B s mixing and B — > irK decays, the final allowed ranges 
for the Z' couplings are presented in Figs. QjandHJ As shown in Fig. (U the range of <\) L S in SI 
is now further restricted with the constraints from B — > irK decays included (i.e., the range 
<frs > —50° is now excluded), while their effect for S2 case is tiny. Our numerical results for Z' 
couplings are summarized in Table. HI 

For evaluating B — > ttK decays, since the Z' mediated effects can occur not only in the 
coefficients of the electro-weak penguin operators but also in the strong penguin ones, we do not 
assume B dd / B^ u = ed/e u = —1/2 (if assumed, the Z' effects in terms of the Wilson coefficients 
of the SM QCD penguin operators will then vanish, as usually adopted in the literature, see for 
example Refs. [131 [HI [15]). So, compared with \B dd \ = \ — B^ u /2\ < 0.1 [H], our fitting result 
in Fig. [2] shows that a larger ranges for B^ u dd are still allowed. Furthermore, as can be seen 
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from Fig. [2] and Table. HJ B dd (B^ u ) is a bit larger than B^ u (B dd ). However, remembering that 
the Z' contribution to the electro-weak penguin coefficient ACV is the combination B^ u — B dd , 
our numerical result always satisfies \B^ U — B dd \ < 0.13, which is still consistent with the result 
31-B^I < 0.3 obtained very recently by Barger et al. 



3 Constraints on Z' couplings from B s — > jj, + fi and B — > 



3.1 Theoretical Framework 

Within the SM, after dropping the negligible charm contributions, the effective Hamiltonian 
for purely leptonic B s — > decay is given as 



H s Jj(B s l+l- 



a 



y/2 2vr sin 2 6 W 



V tb V t * s Y(x t )(sb) v _ A (ll) v _ A + h.c. 



(16) 



where a = -f^ = 1/137, sin 2 Q w = 0.23119 [SB], and the function Y{x t ) is defined as 



O.s 



Y( Xt ) = Y ( Xt ) + -r^Y^Xt) , 



Yo(x t ) 
Y^xt) 



X t r Xt - 4 



3xt 



-x t - 1 (x t 
Ax t + 16x 2 + Ax^ 2x t + %t T ■ 



(17) 



+ 



3(1 -x t ) 2 (l-x f ) 2 
- 10x 2 t -x 3 t -~ 



Li 2 (l - x t ) + 8x t — In 



4x f — 10x 2 — x+ — xf , 2x t — 14x 2 + x] — x\ 



In x t + 



In x t 



(l-x t ) 3 2(1 -x t f 

Within our approximations for the non-universal Z' couplings, the effective Hamiltonian for 
b sl + l~ transition induced by the new Z' boson could be written as 



" (sb) v . A (ll) v . A -^^(sb) v . A (ll) v+A +h.c. (18) 



v tb v t * s 



v tb v t * s 



Then, the full expression for the branching ratio of B s — > l + l is 



x 



a 



2n sin 2 ^ 



Y{x t ) - 2 



v tb v t % 



(19) 
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To clarify the differences between Eq. (TT^|) and the ones in the literature^, a detailed derivation 
of Eq. ffl~9|) is presented in Appendix B. 

The SM effective Hamiltonian for rare b — > sl + l~ decay at scale /i is given by 

H s ef M f (b -+ sl+r) = -^£v tb V t * s [J2 C i^)Qi + Cwi^Qw + C 10A (fx)Q 10A \ + h.c. . (20) 
Here we choose the operator basis given by Refs. [201 EI], in which 

Q 9V = (sb) V -A(Tl)v , Qwa = {sb) v _ A (Jl) A . (21) 

Modifying the Z'-induced effective Hamiltonian Eq. (I18|) to the above form, we find that the Z 1 
effects can be represented as some modifications of the Wilson coefficient of the corresponding 
operators. To this end, the initial conditions for the coefficients Cgy and Ciq A at the matching 
scale \i = Mw are given as 

Cqv,wa(M w ) = C gvloA (Mw) + AC' gVA0A (Mw) , (22) 
&C*{M W ) = -2-2^(B£ + B*) , (23) 

AC? oA (M w ) = 2 -^L(B£-itf), (24) 

with Cgyf WA and AC' 9V10A (Mw) denoting the SM and NP parts respectively. The RG running 
of these Wilson coefficients has been detailed in Refs. [2T)l I2"T| [33]. 

Introducing the normalized dilepton invariant mass s = (pi+ + pi-) 2 /ml, the differential 
decay rate with respect to s for b — > sl + l~ reads 



±J3(b -> sl+l-) a 2 | 




I 2 (1-s) 2 


B{b -> cl-v) 4vr 2 


\Vd,\' 


1 fixXx) 



with 



9+ 2 ~ 9 9+ 2 

L>(s) = (l + 2s)(l + ^-)|C 9 e// | 2 + 4(l + |)(l + ^-)|C 7 e// | 2 

+ [(1 + 2s) + ^-(1 - 4s)] |C 10 | 2 + 12(1 + ±L)C e 7 ff Re(Cf f *) , (26) 

where t = mi/m b , x = m c / m b and E(B — > X c l~v{) = (10.1±0.4)% [3B]. The phase-space factor 
f(x) and the 1-loop QCD correction factor k(x) for B — > X c l~u decay are given respectively 



1 A sign in Eq. (Bl) in Ref. [T3] and Eq. (15) in Ref. [35] mistyped, and the interference terms missed in 
Eq. (Bl) of the first version Ref. [14]. We thank T. Liu and C.W. Chiang for confirmation 
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by 



fix) 



^ 2 + 8x 6 -X 8 -24 X 4 ln X 



1 - 



3tt 



(- 2 -y)(i-x) 2 + 



(27) 
(28) 



The effective coefficient Cg^ is denned as 



C 9 ij(s) + h(x, s)(3C7i + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) - -/i(l, s)(4C 4 + 4C 4 + 3C 5 + C 6 ) 



■i/i(0, s)(C 3 + 3C 4 ) + jj(3C 3 + C 4 + 3C 5 + C 6 ) 



(29) 



where the function rj(s) in the first term represents one gluon corrections to the matrix element 
of Qw, the other terms arise from the insertions of four-quark operators (indicated by the Cj) 
to the one-loop matrix element of Qgy [2H EH [39] 



r){s) 



1 + — 7T 2 Li 2 (s) n s ln(l - s) : r ln(l - s) 

Ti 1 9 3 V ; 3 V ; 3(1 + 2s) V ; 

2s(l + s)(l-2s), A 5 + 9s -6s 2 , 
■his + — — . — ] , 



3(1 - s) 2 (l + 2s 
m b 



6(1 - s)(l + 2s)- 
4 



(30) 



— In In y H 1 — x 

9 /i 9 * 27 9 

2 i I In I ^^- + ] for x = ^ < 1 

— (2 + x)|l-xM l v^F-ii 

^ 2 arctan J— r for x 



^>1 



/i(0,s) 



8 8, m 6 4. . 4. 

In In s H — 27T . 

27 9 /I 9 9 



(31) 



(32) 



Besides these well defined short distance contributions, there are long distance corrections 
related to the cc intermediate states. Phenomenologically they are estimated with Briet-Wigner 
approximation [4"0"| |4"T} |4"2"] which results in a modification to by 



^(3C 1 + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) £ r(K ^ /+r)m ^ 



tt: 



m v n ~ sm 2 b - im Vn T Vn 



(33) 



Usually the factor k ~ 2.3 (more precisely 2.3 times an arbitrary strong phascljH] ) is intro- 
duced phenomenologically to include the known factorizable and the unknown nonfactorizable 
contributions to account for the present experimental data on B — > ^/(nS)X s decays. This 



2 we thank the referee for bring this point to us 
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approximation, however, will cause a double counting of partonic and hadronic degrees of 
freedom[33]. To avoid the double counting, it has been suggested that the long-distance effects 
could be estimated by means of experimental data on R c (s) = a(e~e + —> cc) /a(e~e + — > 
using a dispersion relation |43j(KS approach). In KS approach, only factorizable effect i.e. the 
cc in color singlet, could be estimated with R c (s). It still needs the phenomenological enhance- 
ment factor k to model possible nonfactorizable effects to match the aforementioned large rate 
of B — y ty(^')X s . However, as discussed in detail in Ref. [H], unlike the cc contribution to 
the e~e + — y hadrons cross section where the imaginary part of a current-current correlator 
is integrated over phase space, the huge charm-resonance contributions to B — y X s £ + £~ are 
related to a drastic failure of quark-hadron duality in the narrow-resonance region for integrat- 
ing the absolute square of the correlator over the phase space. Further detailed discussions 
on these non-perturbative effects could be found in Ref. jB]. In this paper, we concentrate 
on the short-distance effects and use the low- and the high-s data, i.e., away from the ^> and 

peaks, to constrain NP effects, and ignore the resonance effects. It should be noted that, 
unlike \l/ and (large data samples, known structures, etc.), the backgrounds due to higher 
J PC = 1 charmonium resonances in the high-s region may be very hard to be vetoed ex- 
perimentally. Although their effects are expected to be much smaller than the former ones, 
they still cause sizable uncertainties which are hard to be estimated [12] • In a recent study of 
exclusive B —y K£ + £~ decay[l6], it is argued that duality violations from the higher resonances 
in high- s region are at a moderate level which may spoil the precision of theoretical predictions 
for (partially) integrated branching ratios of B — y K£ + £~ at the level of several percentage. In 
view of the large uncertainties included in our numerical analyses, we may expect the effects of 
higher cc resonances would not alter our conclusion much. 

Finally, the normalized forward-backward (FB) asymmetry distribution is defined as 



A FB {s) 



li dz &-I-i dz 
Jo dz lStz + I-i dz 



d 2 r 

dsdz 



dsdz 




U 2 E{s) 



(34) 



with 



E(s) = Re(C e / f C* 10 s + 2C e / f C* w ). 



(35) 
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Figure 3: The allowed regions for the parameters in Case I. 



3.2 Numerical analyses and discussions 

With the relevant theoretical formulas collected in section 13.11 and the input parameters sum- 
marized in the Appendix A, we now proceed to present our numerical analyses and discussions. 
The rare B — > X s fi + fi~ and B s — > fi + fi~ decays involve not only the coupling B^ b , which has 
been severely constrained by B s — B s mixing and B — > irK decays discussed in section El but 
also the unrestricted ji — ji — Z' couplings, B^ and . So, our analyses are further divided 
into the following three cases with different simplifications for our attention, namely 

• Case I: with Bfj arbitrary, while taking = 0; 

• Case II: with B^ arbitrary, while taking B^ = 0; 

• Case III: with both B^ and B^ arbitrary. 

In the following discussions, we quote the fitting results for \B^ b \ and 0f under the constraints 
from CB q i <Pb s (95% prob.) and B — > irK decays as inputs. In our numerical evaluations, we 
don't consider the LD contribution to Cg . In each case, our fitting is performed with the 
experimental data on B(B — > X s fi + fi~), B H (B — > X s fi + fi~) and B L (B — > X s fi + fi~) varying 
randomly within their respective la error bars, while the theoretical uncertainties are obtained 
by varying the input parameters within the regions specified in Appendix A. Moreover, we leave 
B(B S — > fi + fi~), A FB (B — > X s fi + fi~ ~) and Ap'g(B -> X s fi + fi ) as our theoretical prediction, 
which could be tested by more precise measurements in the coming years. 
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Table 5: Numerical results for the parameters B^ and B^ (in unit of xlO 2 ). 



Cases 


Case I 


Case II 


Case III 




SI S2 


SI S2 


SI S2 


B L 
mm 


-2.5 ±2.7 -0.55 ±1.0 




-2.7 ±2.5 -0.59 ±0.93 


B R 

MM 




0.78 ±2.0 0.23 ±0.97 


0.61 ±2.4 0.19 ±0.88 
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Figure 4: The dependence of dB{B -»• X s /j,+ n~)/ds on at s = 0.3, 0.5 and 0.7 with 

|.B^| = 1.09 x 10- 3 (2.20 x 10~ 3 ) and <ft = -72°(-82°) in Case I SI (S2) and Case II SI (S2). 



Case I: with B^ arbitrary, while taking B^ = 0. 

In order to investigate the effects of B^ , we neglect the Z' contributions involving B^ in this 
case. Corresponding to the two solutions SI and S2 for B^ b , we obtain two allowed regions 
for B^ as shown in Fig. |3j and the corresponding numerical results are listed in Table |5j 
Our predictions for B(B — > X s fi + fi~), Afb(B — > X s /i + /i~), including the results at both low 
and high q 2 regions, and B(B S — > are given in Table O Due to the fact that the SM 

predictions, B(B -> X s fi + fx~) = (5.0 ± 0.3) x 10~ 6 , B L (B X s fi + fi~) = (1.8 ± 0.1) x 10~ 6 , 
B H (B X.H+IT) = (0.45 ±0.06) x 10" 6 , and B(B S -»■ n+ff) = (3.1 ±0.2) x 10~ 9 , agree quite 
well with the experimental measurements as given by Eqs. (T4|), (jSJ), ©, and (J?]), the parameter 
space with B^ ~ is still allowed. However, as shown latter, nonzero B^ may have significant 
impacts on the other observables. 

With the central values of theoretical input parameters, Fig. |4] (a) shows the dependence 
of dB(B — > X s fi + fi~)/ds on B^ at different s, from which we can see that the minimal value 



15 




of the differential decay rate appears at ~ —0.03 (—0.005) in SI (S2). Furthermore, the 
dependence of dB(B — > X s fi + fi~)/ds and dA FB (B — > X s fi + fi~)/ds on s, with different values 
of Bp , is shown in Fig. From Fig. (a), one can find that dB{B — > X s fi + fi~) / 'ds is reduced 
at B^ ~ -0.03 but enhanced at B^ ~ -0.06 and 0.02 in SI. In S2, with B^ = -0.02, 
the Z' contributions induced by B^ could enhance dB(B — > X s fi + fi~)/ds. However, as shown 
in Fig. (c), it nearly can't reduce B(B — > X s fi + [i~). So, if the future refined experimental 
measurement on B(B — > X s fi + fi~) is significantly smaller than the SM prediction, S2 will be 
excluded first. For both SI and S2 cases, the Z'-induced effects on dApB^B — > X s fi + fi~)/ds 
are tiny as shown in Figs. [5] (b) and (d). 

In our fitting, we find that the constraints on Bzjr are dominated by B^ H \B X 
while the constraint from B{B S — > is very weak due to the fact that there exits only upper 

bound at the moment. At the quark level, since both B — > X s fi + fi~ and B s — > [i + [i~ involve 
the same b — > sfi + fi~ transition, it is interesting to see the Z' contributions to B s — > 
within the parameter spaces constrained by B — > X s fi + [i~ . 
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Figure 6: The dependence of B(B S —?■ 
and -60° (-75°) with \B^ b \ = 1.09 x 1(T 3 (2.20 x 10~ 3 ) in SI (S2). 



-80° (-85°), -70° (-80 c 



Table 6: Numerical results for B(B -» /i+ fi~)(xl0~ 9 ), B L ' H (B -> X s /i+/i~)(xl0~ 7 ), and 
A FB ( J B^X s /i + /i-)(xl0^ 2 ). 



Cases 



Exp. 



SM 



Case I 
SI S2 



Case II 
SI S2 



Case III 
SI S2 



< 47 3.1 ±0.5 3.0 ±1.1 3.3 ±0.8 3.1 ± 1.4 3.5 ±1.1 3.2 ±1.5 3.5 ±1.1 



43; 



16 



-12 
I- 13 

-4.8 
+5.2 



50 ±7 46 ± 10 49 ± 7 50 ± 7 49 ± 7 46 ±11 49 ± 7 
18 ±3.2 17 ±4.3 18 ±3.2 18 ±3.1 18 ±3.1 16 ± 4.5 18 ± 3.2 



B(B S -> fi+fi-) 
B(B^X S ^+^-) 
B L (B -> X fl /i+/i-) 

B H (B X s n + fi-) 4.4 ±1.2 4.5 ±0.6 4.4 ±1.2 4.4 ±1.2 4.5 ±1.1 4.4 ±1.1 4.4 ±1.2 4.4 ±1.1 

-> X a fi + fj,-) 28 ±0.1 29 ± 2 29 ± 1.2 22 ± 6.6 22 ± 5.6 18 ± 13 23 ± 6.7 

Afb(B ~~ ^ X s fi + [i~) 0.5 ±0.3 1.1 ±1.2 1.0 ±0.8 0.2 ± 0.7 0.2 ± 0.5 0.2 ± 2.0 0.4 ± 0.9 

A% B (B -> Xsfi+fi-) 17 ±1.3 17 ±1.4 17 ±1.4 14 ± 4.4 15 ± 3.9 11 ± 7.3 14 ± 4.5 



Fig. |6] shows the dependence of B(B S — > ) on — Bf? at different <p^. In SI, with 
= 0, we find that B(B S — > is easier to be reduced with a smaller Numerically, 

with 0f ~ -65° and B^ ~ -0.05, we get £(£ s -»■ = 2.5 x 10~ 9 , which is 19% 

smaller than the SM prediction 3.1 x 10~ 9 . However, since \<f>s(S2)\ > \$(S1)\, B(B S -> 
is sensitive to seizable \B^ — Bj} ' \ and could be enhanced for most parameter space of 

<P L S - \ B w ~ B w\- For S1 ( S2 )> with <ft = - 80 ° (- 86 °)> B w = -0.06(-0.02), we find that the 
Z' contributions in case I could enhance B(B S — > ji + by about 18% (12%) compared with 
the SM prediction. 
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Figure 7: The allowed regions for the parameters and B R 
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Case II: with arbitrary, while taking B^ = 0. 

Taking B^ = 0, we are going to evaluate the Z' effects induced by B^ . The allowed regions 
of the Z' parameters are shown in Fig. UJ From Fig. |4] (b), which shows the dependence of 
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B(B — > X s fi + fi~) on Bf[ , one would observe that the minimal B(B — > X s fi + fi~) corresponds 
to the point Bp = 0. So, the Z' contributions induced by Bp are nearly can't reduce B(B — > 
X s fi + fi~), which is confirmed by Figs. [8] (a) and (c). 

In Case I, from Figs. (b) and (d), we find dAps{B — > X s fi + fi~)/ds is not very sensitive 
to B^ . However, comparing Figs. [8] (b) and (d) with Figs. [5] (b) and (d), we find dAps^B — > 
X s n + n~) / ds is very sensitive to the Z' contributions induced by Bp . Since the Z' contributions 
to E(s) are smaller than that to D(s), Aps{B — > X s fi + fi~) can be reduced easily rather than 
enhanced. With Bp = 0.03 (0.015) and the central values of the other theoretical parameters, 
Afb(B — > X s fi + fi~) could be reduced by a factor 12% (9%) in SI (S2) compared to the SM 
prediction. At high/low q 2 region, it could be reduced by 12% (13%)/17% (55%). 

Taking Bj^ = 0, Fig. M shows the dependance of B{B S -> fj+fT) on B^. At B^ = 0.03 
with (f)g — = 65°, we get the small B(B S — > /i + /i~) ~ 2.5 x 10~ 9 , which is 18% smaller than the 
SM prediction. With Bj^ = 0.03, (j) L s = -80°, B{B S ->■ /x + /x - ) is enhanced by about 10%. 

Case III: with both B^ and B^ arbitrary. 

More generally, we give up any assumptions for the Z 1 couplings B^ and B^. Because of 
the interference effect between B^ and B^ , the allowed regions for these two parameters are 
now larger than the ones in Case I and Case II, which are shown in Fig. [91 From the figure, 
one can find the correlations between the coupling parameters. Since B^ b and its phase have 
been constrained by B s — B s mixing, B^ and B^ are found to be small. It is interesting to 
note that a model independent constraint on C^ p and it phase has been performed in Ref. [17] 
with Afb(B — > X s ££). Combing our constraints on B^ b , B^ and B^, and re-scaling the 
combination by V ts , one can find our constraints are in good agreement with the magnitude of 
C± P in Ref.JUJ, but with a much stronger constraint on its phase. 

The dependence of dB(B — > X s fi + fi~)/ds and dApB^B — > X s fi + fi~)/ds on B^ and B^ has 
been discussed separately in the last two cases. So, including all of the theoretical uncertainties, 
we just present the dilepton invariant mass spectrum and the differential normalized forward- 
backward asymmetry in Fig. [TUJ We find that in SI, as shown in Fig. [10] (a), BiB — > X s fi + fi~) 
could be either enhanced or reduced by Z 1 contributions. Since the Z' contributions induced 
by Bj^ can hardly reduce B(B — > X s n + n~) as discussed in Case II, the Z' effects for reducing 
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Figure 9: The allowed regions for the parameters and B 



B(B — > X s fi + fi ) are dominated by as discussed in Case I, while the enhancement is due 
to both Bff„ and B* 

However, in S2, as shown in Fig. [TO] (c), B(B — > X s fi + fi~) can hardly be reduced, which 
confirms our analysis in Case. I and II. In both SI and S2, as shown in Figs. [10] (b) and (d), 
A FB (B — > X s fi + fi~) can be easily reduced but hardly be enhanced by Bp . 

With both B^ and B^ included, the branching ratio for B s — > is affected by B^ — 

B^ . Numerically, with Bp — Bp = —0.05, </>f = —65° and the central value of the other 
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cIAfb(B — > X s fj, + fi )/ds. The other captions are the same as in Fig. [5] 



theoretical inputs, B(B S — > is reduced by about 19%, which is the same as in Case I. 

With = —0.05, = 0.03, 0f = —80° and the central value of the other theoretical inputs, 
we find that B(B S — > fi + n~) = 4.5 x 10~ 9 , which is about 46% larger than the SM prediction. 
So, if the coming measurements at LHCb and super B factories present B(B S — > ~ 10~ 8 , 

the family non-universal Z' model will suffer a serious challenge. 

4 Conclusion 

In conclusion, motivated by the observed B s — B s mixing phase anomaly and the so-called "ttK 
puzzle" , we have studied a family non-universal Z' model to pursue possible solutions. With 
the constrained b — s — Z' coupling by B s — B s mixing and B — > ttK decays, we focus on the 
Z' effects on the rare B — > X s fi + fi~ (including both the high and the low q 2 regions) and the 
purely leptonic B s — > fi + fi~ decays, both of which are also induced by FCNC b — > s transitions. 
Our main conclusions are summarized as: 
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• B s — B s mixing anomaly and u irK puzzle" could be moderated simultaneously within 
such a family non- universal Z' model. Corresponding to the two fitting results SI and 
S2 by UTfit collaboration, a new weak phase 0f ~ —72° and —82° are crucial to resolve 
these two problems. 

• Similar to the hierarchy of the CKM elements |^Vt b /y t *^ fe | ~ 0.2, we find \B db /B sb \ ~ 
(^(lCr 1 ) (< 0.2). So, such a hierarchy should be hold within the model. Our results also 
imply the relations B^ u < B L dd and B* u > B* d . 

• Combing B^ b restricted by B s — B s mixing and B — > ttK decays, and B^jf 1 by B — > 
X s fi + fi~, we find B^ R ~ C(10~ 2 ). For observable 8(B — > X s fi + fi~), the reduction 
effects is dominated by the Z' contributions induced by Bj^ in SI. And, both the Z' 
contributions induced by and Bff are helpful to enhance it. The forward-backward 
symmetry Afb(B — > X s fi + fi~) is sensitive to the Z' contributions induced by B^ but 
dull to the one induced by . 

• With the strictly constrained Z' couplings by B s —B s mixing, B — > ttK and B — >■ X s /j, + fi~ , 
comparing with the SM prediction, we find B(B S — > could be reduced/enhanced 
about 19%/46% by Z' contributions at most. The minimal value of B(B S — > 
appears at the point B^ — Bjf ~ —0.05 with the minimal new weak phase <j)^ ~ —65°. 

The refined measurements for the (semi-)leptonic B^ decay in the upcoming LHCb and 
super B factory will provide a fertile testing ground for the SM and possible NP. Our analysis 
about the Z' effects on the observables B^ H,L \B — > A^ L \B — > X s /i + /i~) and 

B(B S — > are helpful to confirm or refute the family non- universal Z' model. 
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Appendix A: Theoretical input parameters 

For the CKM matrix elements, we adopt the fitting results from UTfit collaboration [3J 

p = 0.154 ± 0.022 (0.177 ± 0.044), rj = 0.342 ± 0.014 (0.360 ± 0.031), 
\Vtd/V t s\ = 0.209 ± 0.0075 (0.206 ± 0.012), 

\V cb \ = (4.13 ± 0.05) x 10~ 2 ((4.12 ± 0.05) x 10" 2 ), (36) 

with p = p (1 — 4r) and fj — rj (1 — 4r). The values given in the brackets are the CKM parameters 
in presence of generic NP, and used in our calculation when the Z' contributions are included. 

As for the quark masses, there are two different classes appearing in our calculation. One 
type is the current quark mass which is scale dependent. Here we take 

m s{p) nn A \ n A VJm — /r>r^„\r\ o-7 i CA/r^r n7il — /— \ i ot+0.07 , 



27.4 ± 0.4 g9], m s (2GeV) = 87±6MeV g9], m c (m c ) = 1.27^{ GeV [36] 



m q (p) 

m b (m b ) = 4.20lg;J? GeV [SB] , m t {m t ) = 164.8 ± 1.2 GeV [3S] , (37) 

where m q {p) = (fn u + m^)(/i)/2, and the difference between u and d quark is not distinguished. 
The other one is the pole quark mass. In this paper, we take [361 150] 

m u = nid = m s = 0, m c = 1.6ljT 01 | GeV, 

m b = 4.7918:51 GeV, m t = 172 A ± 1.22 GeV. (38) 
As for the B-meson lifetimes and decay constants, we take [361 EI] 

r Bu = 1.638 ps , r Bi = 1.530 ps , (39) 



f B , = (190 ± 13) MeV , yJB B J Bi = (216 ± 15) MeV 



f Bs = (231 ± 15) MeV , yjB B J B , = (266 ± 18) MeV . (40) 
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Appendix B: Derivation for the Eq. (I19l) 



From the effective Hamiltonian for B s — >■ l + l decay given by Eqs. ( 1T61) and (fl8j) . the amplitude 
for the B s — > l + l~ decay can be written as 



A = A SM + A Z <, 



Asm 
A z , 



G } 



a 



Gp 

71 



(41) 

Y(x t )(^ + fi-\sY(l -75)6 8)^(1 -iMBs) , (42) 



v tb v* 



^ (MVl^l " 7s)& ® j*7„(l - 7 5 )^> 

V tbV ts 



- ' 2 ^f^ (/^Vl^(l " 7s)& ® ^7m(1 + 75>|£.) 



(43) 



a 



After parameterizing the hadron parts, the above equations could be rewritten as 

-F(x t )/i7 M (l -75K (44) 

(45) 



Asm = ^V tb VZif B Pl 



G, 



Az> = ^v tb v; s if Bs p^ 



2% sin 2 6 W 
-2B L M L 



v tb v t * s 



B^ 

1M £7m(1 - 75)/^ + T/ I^ ^mI 1 + 75)^ 



where we have defined (0|s7 M (l — 7 5 )6|i? s ) = —if Ba P^ s , with / Bs being the _B s -meson decay 
constant. For simplicity, we introduce 



Ai = (J,P Bs {l -75)^, 
B = ^V tb V*if Bs , 



fs 

-2B^B L „ 



/</' 



A 2 

c 

Do = 



^B.(l 

n 



2tt sin 2 



75 V, 

-n^t) 



—2BhB 



w 



sb fifi 



v tb v t * s 



(46) 



Then the total decay amplitude can be written as 



A = B[(C + D l )A 1 + D 2 A 2 ], (47) 
\A\ 2 = \B\ 2 [\C + D^ 2 ^ 2 + \D 2 \ 2 \A 2 \ 2 + (C + D 1 )*D 2 A* 1 A 2 + (C + D 1 )D* 2 A 1 A* 2 ] .(48) 



It is easy to get 



IAI 2 
A\A 2 



\A,\ 2 



5 2 2 



A X A\ = -8m 2 m| s 



(49) 
(50) 
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So, 



\A\ 



B\ 8m m B \(C + Di) — D 2 \ 



^\V tb V;/8mlml a 



a 



27rsin 8\y 



Y(x t ) - 



fifi) 



V tb V t * s 



Finally, with \P C \ — f y — 4m^, we get 
B(B s -)/iV) 



iirm 



B s 

r 2 

u F r2 2 

rB s ^f Ba m^m Ba 



m 



x 



ft 



27rsin 9w 



Y(x t ) - 2 



B ' 



V tb V t * s 



(51) 



(52) 
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